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1. Introduction 

In the last years lot of efforts have been devoted to extract interquark 
potential from lattice gauge theories (LGT's) looking at the expectation 
values of Wilson loops or (correlations of) Polyakov lines in Montecarlo sim- 
ulations. Besides the important issue of obtaining reliable values of physical 
observables like the string tension, these simulations also allow to study the 
physical nature of the potential. In particular, more and more convincing 
evidences have been reported in last years of the so called "string picture" 
of the interquark potential: quark and antiquark linked together by a thin 
fluctuating flux tube 

Roughly speaking, there are two independent ways to check this picture. 

The first one is to try to observe directly the flux tube, for instance as 
a stable semiclassical configuration, surviving cooling |^] or looking at the 
space distribution of colour fiux |3| . 

The second one, followed in this paper, is to look at the finite size eff'ects 
due to quantum string fluctuations, which, in finite geometries, give signi- 
ficative, measurable contributions to the interquark potential. This second 
approach traces back to the work of Liischer, Symanzik and Weisz Q and 
has interesting connections with the conformal field theory (CFT) approach 
of two-dimensional models developed in the last years|p. 

This approach has particularly relevant consequences if one studies Po- 
lya_kov lines in LGT's at finite temperature, because in this case the finite 
size corrections (which we shall call from now on "string contributions") 
have a very specific behaviour, due to the choice of boundary conditions. 

The existence of the string contributions and their explicit asymptotic 
form can be deduced from a purely theoretical point of view as consis- 
tency conditions of non-perturbative descriptions of the infrared behaviour 
of gauge theories. 

These corrections seem necessary to fit the data of Montecarlo simula- 
tions [P, 0]. In order to give further evidence of this fact, we shall compare 
our theoretical predictions with a new set of high precision data on the 
Ising gauge model at finite (but low) temperature. In particular, taking into 
account the non trivial topology of the finite temperature model, we shall 
be able to estimate the conformal anomaly c of the CFT which describes 
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the long distance behaviour of the transverse displacements of the string : 
it turns out to agree with the value c = 1 suggested by the naive bosonic 
string picture. 

This paper is organized as follows: after a general introduction to the fi- 
nite temperature gauge models (sect. 2), in sect. 3 we motivate our approach 
and show the general features of contributions coming from the quantiim 
fluctuations of the flux tube. Section 4 is then devoted to the explicit de- 
scription of the finite temperature string corrections . In sect. 5 we describe 
our simulation and discuss the results. 

2. Finite temperature gauge theories: general setting and no- 
tations 

The partition function of a gauge theory with gauge group G regularized 
on a lattice is 

Z = j^dUi{x,t)e^^{-(3Y,ReTr{l-Up)} , (1) 
•> p 

where Ui{x,t) G G is the link variable of position {x,t) and direction I € 
{xi, ..Xd, t} and Up is the product of the links around the plaquette p. 

Let us call Nt (Ng) the lattice size in the time (space) direction (we 
assume for simplicity Ns to be the same for all the space directions). A 
{Ns fNt lattice can then be interpreted as representing a system with a 
finite volume V = (Nsa)'^ and a finite temperature T = l/L = 1/Nta where 
a is the lattice spacing. Lattice simulations with non-zero temperature are 
obtained imposing periodic boundary conditions in the time direction. 

The order parameter of the finite temperature deconfinement transition 
is the Polyakov line, i.e. the trace of the ordered product of all time links with 
the same space coordinates and it is closed owing to the periodic boundary 
conditions in the time direction: 

Aft 

P{x)=TrX{Ut{S,z) . (2) 

z=l 

The vacuum expectation value of the Polyakov line is zero in the confining 
phase and acquires a non-zero expectation value in the deconfined phase. 
The value PdT) of this deconfinement transition is a function of the tem- 
perature, and defines a new physical observable Tc. 
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The interquark potential can be extracted by looking at the correlations 
of Polyakov lines in the confined phase. The correlation of two lines P{x) 
at a distance R and at a temperature T = 1/L = 1/Nta is given by 

{P{x)P^x + R)) = e-^^^'^^ , (3) 

where the free energy F(R, L) is expected to be described, as a first approx- 
imation, by the so called "area law" : 

F{R,L)r. F{R,LU = aLR + c{L) , (4) 

where c(L) is a constant depending only on L. 

The observable @ is similar to the expectation value of an ordinary 
Wilson loop except for the boundary conditions, which are in this case fixed 
in the space directions and periodic in the time direction. The resulting 
geometry is that of a cylinder, which is topologically different from the 
rectangular geometry of the Wilson loop. 

3. The role of the quantum fluctuations of the flux tube 

According to an old, long standing conjecture [^, the area term is only 
the dominant (classical) contribution of the effective string which should 
describe the infrared behaviour of the gauge theory (this explains the sub- 
script "cZ" in eq. (I)). The effective, non critical strings are (as it is well 
known) intractable in more than one dimension. Mainly as a consequence of 
that, it has become a common habit in LGT's to ignore the quantum correc- 
tions implied by the effective string hypothesis, only retaining the classical 
(area) contribution in fitting the data of Montecarlo simulations. We shall 
see below the range of validity and the consequences of this "zero order 
approximation" . 

A way out of this problem is to realize that, for large quark separations, 
we do not need to describe the full complexity of the effective string theory, 
but only the degrees of freedom which survive at large distance. 

Indeed the notorious troubles of the non critical strings ( the breaking of 
the Lorentz invariance or the coupling to the Liouville modes) asymptotically 
disappear at large distances In such a region the transverse degrees of 
freedom of the string are described by the free, massless bosonic fields of an 
ordinary two-dimensional CFT. 
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The corresponding partition function Z{R, L) will be the first order 
quantum contribution of the effective string to the Polyakov line correla- 
tion (or Wilson loop) expectation value. Consequently, taking into account 
the corresponding free energy F{R,L)g = - log Z{R,L), we modify eq.(|) 
as follows: 



F{R,L)^ F{R,L)^i+F{R, L)g = aLR + c{L) -log Z{R,L) . (5) 

Due to the recent progress in CFT's |^] we can study the behaviour 
of F{R, L)q as a function of R and L in a general way. Indeed any two 
dimensional CFT is completely described once the conformal anomaly c, the 
operator content hi and the operator product algebra (or the fusion algebra 
which equivalently encodes all the fusing properties of the CFT) are given; 
then it is easy to show that F{R, L)q only depends on the adimensional ratio 
z = 2R/L. It is possible to give asymptotic expressions for F[R, L)q in the 
z ^ 1 and z ^ 1 regimes: 
^»1 PI 

where d is the space time dimensions of the gauge theory and c — c — 24/i^ 
is the effective conformal anomaly [0]; /imm is the lowest conformal weight 
of the physical states propagating along the cylinder. In the case of unitary 
CFT's hmin = (unless special boundary conditions are chosen) and c 
coincides with the conformal anomaly c; 



F(i?,L),~ -(d-2) c— , (6) 



Z < 1 11 



F{R,L\^-{d-2)c^ , (7) 

where c = c — 24/i„^^ and /i^,/? is the lowest conformal weight compatible 
with the boundary conditions a and (3 at the two open ends of the cylinder. 

If the CFT is exactly solvable, namely if the whole operator content is 
known, one can explicitly write the free energy for all values of z, which will 
smoothly interpolate between the two asymptotic behaviours. 

In the Wilson loop case, since all the boundaries of the loop are open 
(rectangular geometry), we only have the behaviour given by eq.(^), with a 
partition function F{R, L)q symmetric under the exchange R L. 

An important role in this construction is played by the modular trans- 
formations. All the partition functions can be written as power expansions 
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in q = exp(27rir), with t = iz for Polyakov line correlations or r = for 
Wilson loops. Modular transformations allow to extend these expansions in 
the whole r plane. In particular we will be interested, in the following, in 
the r —^/t transformation. 

In the Wilson case, this transformation is a symmetry, because it ex- 
changes R and L. We can use this symmetry by choosing for instance 
L > R, and r = iL/R. Here L plays the role of time and the interquark 
potential V{R) we want to extract from the data is defined in the limit: 
V{R) = limL-,ooF{R,L)/L. 

In the Polyakov case the situation is different: L and R have a different 
meaning and the modular transformation r — I/t allows us to move from 
the region in which 2R > L to that in which 2R < L. What is new is 
that, due to the modular transformation, in these two regions the string 
corrections have, as we have seen above, different functional forms. While 
in the region in which 2R < L the dominant contribution is, like in the 
Wilson case, of the type 1/R, in the region in which 2R > L the dominant 
contribution is proportional to R, and acts as a finite size correction of the 
string tension. This behaviour will play a major role in the following. 

Let us now look at the other side of the problem: the data coming from 
computer simulation. Indeed, in the last years, as a consequence of the 
improvement in precision and size of the Montecarlo simulations of LGT's, 
serious problems arose in the interpretation of data. These problems can be 
summarized into three items: 

a) scaling deviations, 

b) disagreement between data and the "area, perimeter and constant" 
law, 

c) disagreement between Wilson loop and Polyakov line correlation esti- 
mates of the string tension. 

Let us briefly discuss these items. 

a) It is well known that the most safe way one has to separate physically 
meaningful quantities from lattice artifacts is to test their correct scal- 
ing behaviour near the continuum limit, and this is the only ultimate 
check of the whole Montecarlo procedure. Indeed, one of the earliest 
success of LGT's was the identification of such a scaling behaviour in 
the string tension p^ ]. 

However the improvement in the precision and size of lattice simula- 
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tions led to the observation of significative deviations from asymptotic 
scaling |13] of the string tension extracted from Wilson loops. In ||6|, 0] 



we have shown that taking into account string corrections a good scal- 
ing behaviour could be obtained for the string tension extracted from 
Wilson loops both in AD SU{2) and 5C/(3) models |] and in 3D Zsand 
Z5 models 0. 

We shall show in this paper a similar scaling behaviour for the string 
tension a extracted from Polyakov line correlations at finite tempera- 
ture. What is peculiar of the Polyakov case is that the string correc- 
tions for z > 1 act simply as a renormalization of a (see eq.@); as a 
consequence, in this region the scaling behaviour of a can be obtained 
even if one neglects string corrections . 

Owing to the fact that in the actual simulations the size L is always 
rather small, it is enough to cut off few short range (i.e. small R) data 
to reach z > 1 and hence scaling. 

This probably explains why in the last years Polyakov line correlations 
have been so popular in estimating the string tension and have almost 
substituted Wilson loops. 

b) The seeming absence of scaling of the string tension extracted from the 
Wilson loops or Polyakov line correlations at short distance is mainly 
due to a sizable violation of the expected area law of the Montecarlo 
data. As a matter of fact, if one tries to fit also the short range data 
with the area law, the resulting is unacceptably high. Moreover, 
setting a lower cut-off on the fitted data, a definite systematic trend 
can be seen in the resulting string tension, which decreases as the the 
lower bound increases (see fig-l). Let us stress that this behaviour is 
impressively universal: it can be seen in the string tension extracted 
both from Wilson loops or Polyakov line correlations, either in non 
abelian SU{2) and SU{3) gauge theories in four dimensions or in 
the Z2 model in three dimensions Q , . Most probably, this systematic 
trend is the main reason of the fact that, as simulations were performed 
on lattices of larger and larger size, and bigger Wilson loops could be 
reached, the published values of, say, the pure QCD string tension 
became smaller: from the 140 of ref. |l2| in 1980 to the ^ ~ 80 

estimate of ref.s [||, |, |l|, |l6|. 

If one goes on in eliminating short range data, acceptable confidence 
levels are eventually reached when all the data below a suitable thresh- 
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old Ra are ignored. 

Unfortunately this agreement with the area law is only apparent; in- 
deed if we increase Ra further, we still observe a decreasing string 
tension, a trend which is hidden by the error bars (linearly growing 
with the distance) thus allowing apparently good (see fig.l). How- 
ever the values of a extracted at these values of R^ do not scale as 
a universal function of /? as they should do in the asymptotic scaling 
region. 

Notice that this behaviour of o" as a function of Ru is not a lattice ar- 
tifact, since it scales, as we shall see below. In particular, we shall 
account for this behaviour by considering the contributions of the 
quantum fluctuations of the string. 

Let us conclude this item with a remark on the physical scales in- 
volved in these quantum fluctuations. In principle, one could ignore 
all these complications simply by looking at very large Wilson loops. 
This corresponds, in the continuum limit, to study quark- antiquark 
pairs with large separations. Taking as an example SU(3) in 4L> and 
physical unities, we can ignore string corrections for distances greater 
than ~ 1. fm., but they are definitely relevant at distances (~ 0.2 fm.) 
typical of heavy quarkonia (cc, bb). 

c) In this paper we study the Ising gauge model at finite, but low, tem- 
perature, where all the effects of the deconfinement transition can be 
neglected]^. In this region the string tension measured with Wilson 
loops should be the same of the zero temperature case. In fact, ex- 
tracting a from Wilson loops orthogonal to the time direction, we 
found that its value, (compared with the zero temperature results of 
ref. 00) is unaffected, within statistical errors, by the compactification 
of the time direction. On the contrary, if we extract the string tension 
from Polyakov line correlations, a different value is obtained, unless 

^ These effects of extreme interest for the effective string picture, since the quantum 
string fluctuations are more relevant near the deconfinement point. However one should 
take into account higher order contributions. We have chosen to avoid these comphca- 
tions in order to make more definite the comparison between Montecarlo data and our 
predictions. 

^ActuaUy, owing to the finite size (32^) of the lattice, also the data of ref. should 
be considered, roughly speaking, as non-zero temperature estimates of the string tension. 
To be precise in p| the range of temperatures was from 0.35A at /3 = 0.740 to 0.61A at 
13 = 0.7525, where the scale A is defined in eq.([l3). 
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string corrections are taken into account. Indeed this is a well known 
phenomenon: the string tension extracted from Polyakov line corre- 
lations is systematically lower than that extracted from Wilson loops 
(see for instance The difference between the two observables is 

only in the boundary conditions: these do not affect the area term 
(hence we expect the same value of a), but they change the string cor- 
rections. As we noticed above, in the Polyakov case, for large values 

~' d-2) 
6W 



of the interquark distance we must subtract to a, the term "^^"^ 



which explains the disagreement between the measured values |17]. 
Our analysis on the 3D Ising gauge model agrees with a value c = 1 
for the conformal anomaly. 

The simplest possible choice for the effective string is given by the 
Nambu-Goto action, which is obtained by identifying the flux tube with 
the string. The quantum fluctuations of this string can be described, for 
large quark separations, by a gaussian model with d — 2 free bosonic fields 
that we shall call in the following the free bosonic string. This is the simplest 
possible 2D CFT and can be solved exactly. 

The functional integration over the {d — 2) bosonic fields can be done 
using the C function regularization, and strongly depends on the topology 
of the world-sheet |l8|, 0, |2^, One finds c = c = c = d - 2 for the 



effective conformal anomalies defined in eq.s(P,0). Taking into account this 
contribution, we get a definite improvement in the description of the data 
and the above mentioned discrepancies partially disappear. In particular the 
problem of point c) is completely solved as we shall show in the following. 

Notice that this improvement is quite general and works for gauge theo- 
ries both in three and four dimensions, both with discrete and with contin- 
uous groups. 

The free bosonic string can be considered as a sort of mean field approx- 
imation of the proper effective string description and as such, for instance, 
predicts the same (mean field) critical index = 1/2 near the deconfinement 



transition for all gauge theories (see ref. |22| for a discussion of this point). 

There are three features of the bosonic string which explain in which 
sense it is an approximation of the physical flux tube: 

i) The bosonic string can freely self-overlap, a feature which does not agree 
with our intuition on gauge theory coming both from strong coupling 
expansions ( in a Z2 gauge theory, each plaquette must appear at 
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most once in any strong coupling diagram) and from semi-classical 
descriptions of the flux tube. 



ii) The bosonic string has zero thickness, while the flux tube has a nonzero, 

measurable thickness [^]. 

iii) The bosonic string predicts, for the short distance behaviour of a p^, 



a(i?) = a(oo)^l-(^)' . (8) 

The term that we can obtain using CFT's (see sect. 4) is only the first 
order one in the 1/R expansion of eq.(^). Consequently we expect 
that it underestimates the true contribution and that, once it has been 
subtracted, the decreasing behaviour of (t{R) should be weaker, but 
still present. This is exactly the opposite of what we observe. Both in 
the case of Wilson loops (see fig. 4 of ref.]^]) and in the case of Polyakov 
line correlations (see fig.l of the present paper), once bosonic string 
corrections are taken into account, cr{R) has the wrong shape. 

The free bosonic string is a good description of the effective string at 
large distances and low temperatures: it works for values of R and of the 
temperature such that the finite thickness of the flux tube can be ignored, 
and the self-overlapping probability is negligible. 

A simple way to take into account the presence of a new scale (namely 
the finite thickness of the flux tube) in the problem, is to compactify into a 
circle the bosons describing the transverse displacements . There are several 
ways to look at this compactification: one may either think that the string is 
not any more identified with the flux tube but wraps around it (2^] , or notice 
that a compactified boson can be fermionized and hence one can look at the 
model as a fermionic string description [^,0, 25 1 or, using the connection 



between critical 0(n) models and compactified bosons, one can discuss the 



string corrections in term of self- avoiding walks |26]. 



All these approaches capture some of the physical content of the model, 
and can give some semiclassical and intuitive justifications, but let us stress 
again that the ultimate reason for introducing a compactification radius is 
only to set a new scale in the model. This scale does not represent a new 
degree of freedom, since general consistency requirements fix the compact- 
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ification radius to be exactly 1/4 in conventional units [^] and select the 
operator content of the underlying CFT 26|. 



The compactification on a circle of rational radius leads to a simpler 
description of the CFT in terms of Dirac fermions with a finite number of 
spin structures 1,0, H]. In the following we shall call this model "fermionic 
string" . As expected, the large distance behaviour of the fermionic string is 
the same of that of the free bosonic string. In particular, one can show that 
again c = c = 1, hence both strings predict the same gap between Polyakov 
and Wilson estimates of a. The difference between the two models becomes 
sensible when short distance data (of the order of the flux tube thickness) 
are taken into account: in fact it can be shown that now c = (d — 2)/4. In 
this region the area law corrected with fermionic string contributions gives 
definitely better x^'s, and reliable scaling values of the string tension can be 
extracted (see for instance ref.s [^0]). Finally, a further signature of the 
better behaviour of the fermionic corrections with respect to the bosonic 
ones is given by the shape of the string corrected function (j{R) which, as 
expected, is still a weakly decreasing function of R (see fig.l for the Polyakov 
line correlations and fig. 4 of ref. for Wilson loops). 

In this paper we shall confirm this picture also in the finite temperature 
case. 

4. String corrections in the finite temperature geometry 

The case of the free bosonic string was discussed in ref.s [|^, |l^ and 
with a different approach in ref. ||2^ (while the calculation for the rectangular 
geometry can be found in [^]). The result is, for (i = 3, 

F,(i?,L)=log(7?(r)) ; -ir = ^ , (9) 
where rj denotes the Dedekind eta function: 

oo 

Vir)=q^^l[{l-qn ; 9 = e^^^ , (10) 

n=l 

and R is the distance between the two Polyakov lines. 

We list below for completeness the power expansions in the two asymp- 
totic regions: 
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2R<L 

J oo 

F,{R,L) = -^ + Y.^og{l-e---'^"') , (11) 



2AR 

n=l 



2R> L 



F,(i?,L) = -^ + ilog^ + f:iog(l-e-^-"«/^) . (12) 

n=l 

Following the same line of reasoning, one can also obtain the string 
correction in the case in which the transverse degrees of freedom of the 
string are described by fermionic fields^. We shall list here, for completeness, 
the results for generic values of the compactification radius r, and then 
concentrate on the value r = 1/4. 

The string is described by a Dirac spinor ip{^i,S,2) (0 < < i?, < ^2 < 
L) with components ip±{Ci,(,2) = fp±{^i ±^2), associated to the transverse 
dimension. 

The value of the quantum correction depends on the boundary conditions 
of ip on the two fixed edges of the cylinder (the Polyakov lines) . Defining a 
boundary phase a = through the relation 

^±(6 + 2i?,6) = e^'"*/^V'±(6,6) , (13) 

one has 

Ff,N{R,L) = - log J2 ; -tT = — , (14) 

m=0 ^ 

where t?[^](0|r) is the Jacobi Theta function with characteristic [^] . One 
can write explicitly (0|r) in terms of an infinite product in analogy with 
eq.(0): 

^/^^^ ' =g¥-2i J|(l + g"+°-ie2^^^)(l + g"-"-^e-2^^^) . (15) 



?7(r) 



n=l 



Given the above expressions the series expansions for z < 1 and z > 1 
are, in the case of N=4: 

^This was done for rectangular geometry (Wilson loops) in ref.s ^, Q in the case of 
Majorana fields (Ising model) and for Dirac fermions (c=l conformal field theory) with 
arbitrary boundary conditions. 
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2R<L 



FfA{R,L) 



ttL 



- log[4(l + q + 2q^ + ■■■)] 



(16) 



96i? 



with q 



exp{ 



nL \ 
R ) 



2R> L 



Ff,^{R, L) = - log(l + \l2 + V2q^l^^ + ^/2q^l^ + • • •) (17) 

with q = exp{-^) . 



Notice that while the bosonic and the fermionic strings have quite differ- 
ent behaviours in the region 2R < L due to the different operator content, 
they give exactly the same dominant contribution in the region L < 2R 
since they have the same conformal anomaly c = 1 and both are unitary 
CFT's. In this regime the only relevant difference between the two strings 
is a sub dominant contribution of the type log R: this is present only if the 
boundary conditions of the transverse degrees of freedom are periodic (like 
in the case of the bosonic string) and is forbidden if they are non-trivial (like 
in the fermionic case). 

The above listed contributions are, as a matter of fact, only the first 
order terms of the string corrections to the interquark potential. They are 
expected to be valid in a window far enough from the roughening transition 
so as to allow the string to fluctuate, but not too near to the deconfinement 
transition where the string delocalizes, and the higher order effects due to 
self-overlapping become important. 

5. Comparison w^ith Montecarlo data 

In order to measure the conformal anomaly c we made extensive high 
statistics simulations of a 3D Ising lattice gauge theory, extracting Wilson 
loops and Polyakov line correlations in a wide range of values of the ratio 



where Up denotes the product of Z2 variables associated to the links be- 
longing to the plaquette p. This model is known to have a roughening 



R/L. 



We used the standard action 




(18) 



p 
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transition at Pr = 0.4964 and a zero temperature deconfinement transition 



at (3c{Nt = oo) = 0.7614 [27| which is of second order. The corresponding 
critical index is ~ 0.63 |2^. The scahng region of the finite tempera- 
ture T = l/Nfa deconfinement transition starts approximatively at Nt = 4. 
Data obtained with the Montecarlo renormahzation group approach ||2^ for 
Nt = 4 and Nt = 8 give Pc{Nt = 4) = 0.7315(5) and pc{Nt = 8) = 0.7516(5) 
which correspond to the following value of the critical temperature Tc 

TJ A = T,a{(3,{oo) - (3c{Nt))-'' = 2.3 ±0.1 , (19) 

which is also in agreement with an independent prediction based on the N=2 
supersymmetric behaviour of the critical point pT 



We made six Montecarlo simulations on lattices of sizes {Ns)'^Nt for the 
values of /? and Nt described in tab. I . These values have been chosen in 
order to check both the scaling as a function of /3 and the temperature 
dependence of the physical observables. In fact the first three combinations: 
{(3 = 0.739, Nt = 12), (/3 = 0.742, A^i = 13), {/3 = 0.746, A^* = 15), roughly 
correspond to the same physical temperature T/A ~ 0.92 (see tab. I), while 
the second set: (/3 = 0.739, Nt = 13), (/3 = 0.742, Nt = 14), (/3 = 0.746, Nt = 
16) corresponds to T/A ~ 0.85. These are the lowest temperatures that we 
could reach with our precision. 

The consistency of our data with the zero temperature result tells us that 
they are low enough to be unaffected by the higher order corrections which 
appear near the deconfinement transition. The size in the space directions 
was fixed to be Ng = 48, and periodic boundary conditions were chosen also 
in these directions. A standard heat-bath algorithm was used to update 
links. We measured all the possible Polyakov line correlations, and in the 
two experiments, (/3 = 0.739, iVt = 13) and {(3 = 0.746, A^j = 15),0 we also 
measured Wilson loops (orthogonal to the time direction) of sizes (R,L) in 
the range 2 < R,L < 20. 

The major problem we had to face was the critical slowing down. We 
have already discussed this problem in ref. Q. Comparing the present data 
with those of ref. we can say that the problem is even more severe with 
finite temperature simulations. Moreover, Polyakov line correlations are in 
general more self-correlated than Wilson loops. We kept under control this 
problem separating the measure of two successive Wilson loops and Polyakov 

^We have chosen the two extreme situations, both in f3 and in T. 
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line correlations with 64 sweeps. We also developed a code so as to scatter 
both in space and (Montecarlo) time the measured correlations, trying to 
optimize not only self-correlations but also cross-correlations. We think that 
this procedure could allow a significative improvement in precision for any 
LGT simulation. The typical run was composed by 200.000 sweeps of ther- 
malization and 960.000 sweeps during which 15.000 measures of Polyakov 
line correlations and Wilson loops were taken. 

The resulting expectation values were fitted first with the pure area law 
and then adding the fermionic string quantum contribution (eq.s(p^l7)). 



This correction being fixed, we do not add new free parameters in the fitting 
procedure, hence the two fits (with and without quantum contribution) can 
be directly compared. In fitting the data we had also to take into account the 
periodicity of the lattice: the same correlation between two Polyakov lines 
could also be obtained following topologically non-trivial paths, winding 
around the compactified space directions. We made nonlinear fits summing 
over the nearest winding number^. So the fitting function was 

«,i=i 

{P{x)P{x + R)) = J2 e-^(^»-^'^*) , (20) 
«,i=-i 

with 

= ^{R^ + iN,y + {Ry + jNsY , (21) 

Rx and Ry being the Euclidean components of the vector R and 

F{Rij,Nt) = a RijNt + cost +Ff4Ri,j,Nt) . (22) 

Being Nt fixed, the constant term must be the same for all the correlations. 
Notice that this fit is definitely more constrained than the similar one for 
Wilson loops, since in this case we have only two free parameters, namely 
cost and u, instead of the three of the Wilson loop. 

In fitting the data we used the expansions of Ff^4 up to the sixth order^ 
in q (eq.s(|l6,17)) both in the region z < 1 and in 2; > 1. The Wilson loops 



analysis was done as described in ref. Q. 

^The same problem can be ignored with Wilson loops. 

®As a matter of fact only the dominant term and the second order give significative 
contributions for generic values of z, higher orders become relevant only in the point 
z = 1 where they allow the smooth connection of the two asymptotic behaviours. Since 
this point turns out to be very important in our fits, we decided to keep higher orders 
anyway. 
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Before discussing our results let us make few general comments on the 
data: 

1. A first feature of the Polyakov line correlations is their impressively 
good scaling behaviour (which is shown in fig. 2) even before any cut- 
ting or fitting manipulation. The two sets of data which overlap in 
fig. 2 are taken at different values of /3 and Nt tuned so as to give the 
same physical temperature. The tuning crucially involves the critical 
index = 0.63 and tells us first, that we are in the scaling region and 
second, that, as expected, all the data are controlled only by the two 
physical scales T and a . 

2. As described in sect. 3, we cannot expect to see any string effect in the 
R dependence of Polyakov line correlations at large distances {z > 1) 
(in this region the only signature of the string is in the T dependence 
of a and it can be seen only by comparison with the Wilson loops, 
being the two values of T we explore are too near to give measurable 
effects). On the contrary string effects should be visible for z < 1 as 
deviations from the linear rising behaviour of the potential. This is 
clearly visible in fig. 3 . In this figure the potential V{R), defined as 

ViR) = -^log{P{x)PHx + R)) , (23) 

is plotted against the interquark distance R. The linear rising be- 
haviour is clearly visible for z > 1 (in the example z = 1 corresponds 
to R = 7.5), and similarly evident is the smooth onset of the 1/R 
correction for z < 1. 



Let us now concentrate on the quantitative tests. We made several fits, 
cutting off short range data, as described in sect. 2 and for the Wilson loops 
we cut off short range data following the procedure of ref. |^ . The result of 
these fits is a set of values of the string tension a{(3, Nt, Zmin) (in the case 
of Wilson loops Zmin is substituted by Rmin)-, which depend on /3, on the 
temperature T = l/Nta and on the threshold Zmin {Rmin)- Wilson loops 
data are listed in tab. II together with two sets of data taken from ref.0] for 
comparison. The Polyakov line correlation data are listed in tab. Ill . The 
errors on a have been estimated with an ordinary jackknife procedure, and 
the data are presented taking into account their scaling behaviour in /?: 
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Let us examine first the data coming from Wilson loops. We can see 
from tab. II that: 



If we add the fermionic correction, all data agree among them: for all the 
sets, acceptable confidence levels are obtained already from Rmin = 3, 
and, starting from the same threshold, scaling is fulfilled. 

On the contrary, if we neglect quantum corrections, acceptable confidence 
levels with the area law are reached only at Rmin = 5, and, what is 
more important, scaling is never reached. 

It is quite interesting to see that the data show no dependence at all 
from the finite temperature. This is particularly impressive for the 
results obtained without the string corrections: the differences be- 
tween various /3's due to the scaling violations are much bigger than 
the difference between zero temperature and finite low temperature 
data for nearby values of (3: the data at (/3 = 0.739, T/A = 0.84) 
are much more similar to those at (/3 = 0.740, T = 0) than to 
those at (/3 = 0.746, T/A = 0.92), and similarly for the other pair 
(/? = 0.746, T/A = 0.92) and (/? = 0.745, T = 0). 

The values of a obtained with the string corrections still show a decreasing 
trend in the region Rmin > 3. This behaviour was already noticed and 
discussed in ref. 0: it can be ascribed to higher order contributions 
of the string and can be successfully fitted with a root law The 
important point is that this trend (as a difference with the no-string 
case) becomes asymptotic already in the region Rmin = 5,6 which we 
can reach with our measures. Since the new data are in agreement 
with those at T = 0, we can borrow the result of ref. Q and give 

- (I^Z^ = 3-65(9) (25) 

as a final value for the string tension extracted from Wilson loops. 

Analyzing now the data coming from Polyakov line correlations (tab. Ill) 
we can see that: 

In general, for the string-corrected fits, acceptable confidence levels are 
obtained for z > Zmin = 0.6 , while using the pure area low one finds 
acceptable CL only for z > Zmin = 0.8 . 
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The fits to the scahng law show the same thresholds and give definitely 
good confidence levels for z > Zmin = 1-2, a region in which the string 
correction only acts as a renormalization of a (see eq.(p^). However 
the fact that the fermionic string describes better the data can be still 
stressed observing that, after z = 1.0, the string-corrected data give a 
stable (T, while those obtained with the pure area law keep decreasing. 

Taking as final answer for the string tension the value (and error) extracted 
at Zmin = 1-2, we haveQ 



h = = 3.70(4) (26) 



which is in impressive and remarkable agreement with the value k^^ 
given in eg. (p5|) . 



Let us conclude making the same comparison between Wilson and Po- 
lyakov data, without taking into account string corrections. 

Since in this case Wilson loops data do not scale and Polyakov data 
should depend on Nt, we cannot do this comparison at the level of scaling 
variables and kp as we did above, but we must examine the data at each 
value of (3 and Nt separately. The result is shown in fig. 4 and fig. 5: there 
is clearly a gap between the values of a evaluated with these two different 
observables. 

In order to get a numerical estimate of this gap, we assumed the following 
rule: consider, both for Wilson loops and for Polyakov line correlations, the 
first value of the string tension (as Zmin and Rmin increase) with a good 
confidence level of the fit to the pure area law. As discussed above, this 
means that we must take Rmin = 5 and Zmin = 0.8. With this criterion we 
obtain: 

a{Nt = 13,/? = 0.739 , Rmin = ^)u,iison = 0.0338(7) , 

(j{Nt = 13,/3 = 0.739 , Zmin = ^■^)polyakov = 0.0315(4) , 

with a gap 0.0023(11) ; 

a{Nt = 15,/3 = 0.746 , Rmin = ^)wiison = 0.0214(3) , 

a{Nt = 15, /? = 0.746, Zmin = l-2)polyakov = 0.0195(3) , 

with a gap 0.0019(6) . 

''The absolute value of a does not change within the errors taking any value of 
Zmin > 1.0. Only the error is a function of Zmin and we choose , of course, the z value 
corresponding to the maximal confidence level, which is higher than 70%. 
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Remarkably, both are compatible with a value c = 1 for the effective 
conformal anomaly. Indeed, since these gaps should be due, at least asymp- 

, setting c = 1, we expect a gap of 0.0031 
for Nt = 13 and of 0.0023 for Nt = 15. 

Needless to say that these numbers must be taken with caution, since 
both the Up and are still moving toward their asymptotic values in the 
region that we are measuring. Notwithstanding this, they clearly indicate 
the presence of a gap between the two estimates of a, and this fact can 
simply be understood only in the context of an effective string picture of the 
infrared behaviour of the gauge theories. 



totically, to the term of eq.(17 



Figure Captions 

Fig.l) The string tension cr as a function of Zmin for the fermion, boson 
and no-string pictures at /? = 0.746, Nt = 15. 

Fig. 2) Polyakov line correlations as functions of the adimensional ratio 
z = 2R/Nt. Three samples of data are presented: /3 = 0.746, Nt = 16, 
p = 0.742, Nt = 14 and /? = 0.742, Nt = 13. The first two, which 
overlap in the figure, correspond to the same physical temperature T. 
The third one, with error bars, is plotted for comparison. 

Fig. 3) The inter quark potential V{R) defined function of the 

distance R. The data are taken from the sample at (3 = 0.746, Nt = 
15, error bars are presented only for i? > 8. For R < 8 errors are 
smaller than the plotted symbols. The straight line corresponds to 
the pure area law fit. The dashed line is the fit with the fermionic 
string correction while the dotted line is the fit with the bosonic string 
correction. 

Fig. 4) String tension extracted from the Polyakov line correlations (square) 
and from Wilson loops (diamonds) using a pure area law, without 
string corrections, plotted as functions of the short range cutoff Rmin 
. The sample of data is at /3 = 0.746, A'f = 15. The values of Zmin for 
the Polyakov case are converted using the relation R = z ■ Nt/2. 

Fig. 5) Same as Fig. 4, taking into account the fermionic string corrections. 



18 



Table Captions 



Tab. I Set of measured data: in the third column the corresponding physical 
temperatures are presented. In the last column we list the samples in 
which also Wilson loops were measured. 

Tab. II Values of string tensions extracted from Wilson loops as functions 
of Rm for the two set of data: (3 = 0.739, /? = 0.746 presented in this 
paper and for three sets of data of ref. [^: (3 = 0.740, f3 = 0.745, 
/3 = 0.750. The string tension is presented in scaled unities in the case 
of the pure area law fits (Ila) and taking into account fermionic string 
corrections (lib). In both tables the sixth line contains the result (k) 
of the scaling fit and the last line the corresponding reduced aiid 
confidence levels. 

Tab. Ill Values of string tensions extracted from Polyakov line correlations 
as functions of Zmin- The string tension is presented in scaled unities 
in the case of the pure area law fits (Ilia) and taking into account 
fermionic string corrections (Illb). In both tables the last two lines 
contain the result (k) of the scaling fit and the corresponding reduced 
and confidence levels. 
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Tab. I 



/3 


Nt 


T/A 


Wilson 


0.739 


12 


0.91 




0.739 


13 


0.84 


X 


0.742 


13 


0.92 




0.742 


14 


0.86 




0.746 


15 


0.92 


X 


0.746 


16 


0.87 





l^b. II (a) 



/3 




3 


4 


5 


6 


0.739 


5.03(1) 


4.50(2) 


4.22(4) 


4.06(8) 


4.14(16) 


0.740 


4.90(1) 


4.47(3) 


4.24(5) 


4.05(9) 


4.09(19) 


0.745 


4.97(2) 


4.51(3) 


4.30(5) 


4.17(9) 


3.99(16) 


0.746 


5.14(1) 


4.62(2) 


4.36(3) 


4.12(6) 


4.20(9) 


0.750 


5.30(3) 


4.77(6) 


4.46(8) 


4.32(14) 


4.21(25) 


k 


5.028(5) 


4.55(1) 


4.31(2) 


4.12(4) 


4.14(6) 




35.0 ; 0.0 


10.2 ; 0.0 


3.3 ; 0.0 


0.9 ; 0.5 


0.4 ; 0.8 



Tab. II (b) 
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/3 




3 


4 


5 


6 


0.739 


4.26(1) 


4.01(2) 


3.87(4) 


3.79(8) 


3.92(16) 


0.740 


4.19(1) 


4.01(3) 


3.90(5) 


3.77(9) 


3.72(19) 


0.745 


4.15(2) 


3.98(3) 


3.91(5) 


3.85(9) 


3.73(16) 


0.746 


4.23(1) 


4.04(2) 


3.93(3) 


3.77(6) 


3.89(9) 


0.750 


4.24(3) 


4.07(6) 


3.96(8) 


3.90(14) 


3.84(25) 


k 


4.221(5) 


4.02(1) 


3.91(2) 


3.80(4) 


3.85(6) 




9.7 ; 0.0 


0.9 ; 0.4 


0.5 ; 0.7 


0.3 ; 0.9 


0.3 ; 0.8 
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Tab. Ill (a) 





Nt 


0.2 


0.4 


0.6 


0.8 


1.0 


1.2 


1.4 


0.739 


12 


4.48(3) 


4.04(3) 


3.89(4) 


3.79(5) 


3.71(5) 


3.64(7) 


3.64(12) 


0.739 


13 


4.63(2) 


4.15(3) 


3.93(4) 


3.78(4) 


3.67(6) 


3.57(12) 


3.54(24) 


0.742 


13 


4.58(2) 


4.12(3) 


3.94(4) 


3.83(5) 


3.76(7) 


3.72(10) 


3.69(17) 


0.742 


14 


4.65(1) 


4.16(2) 


3.91(3) 


3.76(3) 


3.64(3) 


3.57(9) 


3.56(14) 


0.746 


15 


4.36(3) 


4.08(3) 


3.83(5) 


3.75(6) 


3.65(8) 


3.62(10) 


3.56(13) 


0.746 


16 


4.52(3) 


4.08(4) 


3.94(5) 


3.82(6) 


3.81(9) 


3.71(13) 


3.79(21) 




k 


4.58(1) 


4.12(1) 


3.91(2) 


3.78(2) 


3.69(3) 


3.63(4) 


3.62(6) 




: C.L. 


10. : 0.0 


:-!.r, : 0.0 


0.7 : ()..") 


0.:-! : 0.9 


0.7 : ()..") 


0.:-! : 0.9 


0.2 : 0.9 



Tab. Ill (b) 



/9 


Nt 


0.2 


0.4 


0.6 


0.8 


1.0 


1.2 


1.4 


0.739 


12 


4.07(2) 


3.89(3) 


3.80(4) 


3.75(5) 


3.71(5) 


3.67(7) 


3.70(12) 


0.739 


13 


4.17(2) 


3.98(3) 


3.84(4) 


3.74(4) 


3.67(6) 


3.59(12) 


3.59(24) 


0.742 


13 


4.14(2) 


3.96(3) 


3.86(3) 


3.81(5) 


3.78(7) 


3.77(10) 


3.77(17) 


0.742 


14 


4.16(1) 


3.97(2) 


3.82(3) 


3.72(3) 


3.63(5) 


3.59(9) 


3.61(14) 


0.746 


15 


4.04(3) 


3.90(3) 


3.77(4) 


3.72(6) 


3.67(8) 


3.67(9) 


3.64(13) 


0.746 


16 


4.15(3) 


3.93(4) 


3.86(5) 


3.79(6) 


3.80(9) 


3.74(13) 


3.81(21) 




k 


-1.11(1) 


:i.95a) 


3.82(2) 


.^7.5^2) 


:?.09(,s) 


:U>7(1) 


:i.(i8(0) 




; C.L. 


4.0 ; 0.0 


1.3 ; 0.2 


0.6 ; 0.6 


0.4 ; 0.8 


0.8 ; 0.5 


0.4 ; 0.8 


0.2 ; 0.9 
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